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Abstract 

Several new results on the multicritical behavior of rectangular matrix models 
are presented. We calculate the free energy in the saddle point approximation, and 
show that at the triple-scaling point, the result is the same as that derived from 
the recursion formulae. In the triple-scaling limit, we obtain the string equation 
and a flow equation for arbitrary multicritical points. Parametric solutions are also 
examined for the limit of almost-square matrix models. This limit is shown to provide 
an explicit matrix model realization of the scaling equations proposed to describe 
open-closed string theory. 
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1 Introduction 



In recent years, the study of random matrix models as a lattice regulator of two-dimensional 
euclidean gravity has shed light on the theory of non-critical strings. The discovery of 
double-scaling limit [1-3] has led to results to all orders in the genus expansion. The study 
of new matrix models may lead to further insights for the continuum theory of random 
surfaces. 

Rectangular N x M matrix models differ from previously studied models in that they 
have two independent large N parameters. In general, one finds that the critical behavior is 
characterized by two-dimensional singularities. There are two special cases in these models, 
described by a double-scaling limit: the vector model limit where M remains finite, and the 
almost-square matrix limit where P = N — M remains finite. These cases were solved in 
ref. [|J, where a surface interpretation was also provided.^ The first limit leads to a phase 
of dense polymers, and the second, to random surfaces. Another class of critical points is 
described by a triple-scaling limit, which corresponds to independently taking both N and 
M to infinity in correlation with the approach of critical matrix couplings^]. 

In this paper, we present new results for the multicritical behavior of rectangular matrix 
models. The remainder of this section introduces our notation. In section ^|, we explicitly 
calculate the free energy for the simplest non-trivial potential in the saddle point approxi- 
mation. This calculation is performed to ensure that the recursion formulae do not lead to 
spurious critical points. As expected, we verify that the triple-scaling calculations correctly 
produce the critical behavior of the free energy. In section |3], we solve the model in the 
triple-scaling limit, and present the string equation for arbitrary multicritical points. The 
partial differential equation arising from the potential independent recursion relations is 
also formulated as a flow equation M. Next, we examine the limit of almost-square matrices 
where P remains finite. We show that in this limit these matrix models provide an explicit 
realization of the open-closed string equations proposed in ref. ||, with P playing the role 
of the open string coupling constant. We briefly discuss the new results in section |5|. 

We now review the techniques developed in ref. |7|] for the solution of rectangular matrix 
models. Given an ensemble of N x M matrices T with complex entries, we want to study 
the partition function 



where V(T^T) = J2p=i dp{T^T) p / (2p). This matrix ensemble is related to the tangent space 
of U(N+M)/(U(N) x U{M)) [§. One can "diagonalize" T with the natural U{N) x U{M) 



where V\ G U(N), V2 G U(M) and Xm = diag (xi, ■ ■ ■ , xm)- Thus the partition function 
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3 The model with N — M was first studied in ref.'s j|, Hi- 
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reduces to 



r°° M r i 
Z oc / exp(-2/SV( yi )) \Vi\ N - M } \{ (vt-Vj) 2 (2) 

J ° i=l l<i<j<M 

where we have set Ui — x\. One observes that the partition function depends independently 
on M and P = N — M. One could try to incorporate the y N ~ M as a logarithmic term in the 
potential, but there is no systematic way to proceed with such a potential. Rather ref. |7j 
developed recursion relations relating measures dfi[l] = dyy 1 exp(— 2(3V(y)) with different 
powers I. This provides a systematic framework to study all phases of these models.^ 
Consider orthogonal polynomials on the positive real line 



The recursion relations stepping in I are 

Pg>(v) = PS +1) (v) + 4 } P^\y) 
yPt x \y) = P^lM + o%iP^\y) ■ 

The above coefficients also enter into the standard recursion relation 

vp2> = p^ + ^I + ^p^ + oXp^ ■ 

Elementary manipulations produce the potential independent relations 

€ - = e x) - *<? e^fc? = m ] ■ ^ 

One may also derive two potential dependent relations from 

2 " + , + > = ldn[i]yV'{y)P^\y)P^{y) 



2(3 

Y/n-? = jdii[l]V\y)[P^\y)P^\y) + e^P^{y)P^ 1 {y)] . (4) 

In the planar limit, (3 — > oo with g = N/(3 and q = M/(3 fixed, we assume QM±j — > 6> and 
VM±j 4> • Then, eq. ([$]) is trivial while eqs. (|J) take the form 

g = ed e u{e,(j>) q = <f ) d^u(e,<f ) ) (5) 

where 

/■27T W\ 

U{6,<f)) = 2 v{e iX + e + <p + e - iX 

Jo 2ir 



J These techniques could also be applied to Penner-like models [nil 
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The critical points are identified as points where the two-dimensional map (g, q) 
is singular. Such points occur at the vanishing of the jacobian determinant 



deg d^g 
deq d^q 







(6) 



Multicritical behavior is produced by demanding higher derivatives of linear combinations 
of g and q also vanish. The free energy may be determined from 
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2 Saddle point approximation 

We begin by calculating the free energy for the simplest non-trivial potential in the saddle 



point approximation using the techniques of ref. [ll]. This calculation provides a test 
to verify that the triple-scaling analysis of the recursion relations correctly generates the 
critical behavior of the matrix model. There are known cases where recursion analysis 
leads to spurious results. For example, using multiple limits (more than two) of the scaling 
functions to describe eigenvalue densities on multiple intervals, generically yields false 



critical points|12[]. Additional constraints producing identical minima for the potential on 



each interval are required. For rectangular matrices in the triple-scaling limit, given that 
different derivations lead to expressions for the free energy which differ by functions of s 
(defined below), one may suspect that the results are also spurious. Below the nonanalytic 
behavior of the free energy identified through the recursion relations is reproduced in the 
saddle point approximation for the simplest non-trivial case. This confirms that the triple- 
scaling analysis of the recursion relations does not lead to spurious results. 

The planar limit of rectangular matrix models with non-vanishing a\ and a2 was first 
solved in ref. [[OJ. For a general potential, the saddle point approximation is discussed 



in ref. 0. We will briefly review the results below. The entire integrand of the partition 
function in eq. (0) may be written exp(— j3E) with 

M 1 M N — M M 

E=2j2V(y i )--j:\n\y i -y j \-——j:iny i . (7) 

i=l P i,J=l P i=l 

Eq. (0) has an interpretation as the energy of M charged particles on a line. The three 
contributions are: an external potential, 2V(yi); a Coulomb interaction between two par- 
ticles; and an electrostatic repulsion away from the origin, with strength P = N — M. At 
equilibrium, the particles are confined to an interval [A, B] with < A < B. 

Following ref. |l]J , the saddle point solution of @ entails the introduction of an eigen- 
value density, p(z) which solves 

dV(y) P M r B dz 



dy 2(3y (3 J a y- z 



M dz 

B JAV — Z 
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One finds that p(z) takes the form 



p(z) = JL u {z)j{B-z){z-A) 



nM 



with u(z) = X)fc=_i where 

L-2-k 



hi 



1 E 
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(9) 



(10) 



After integrating the saddle point free energy becomes 



Fo = (3E = M dyp(y) 2(3V(y)-P\n\y\-2M\n\y-z \ + /3MV(z ) 



PM 

— lnz 
(11) 

where z Q is a constant of integration, which does not affect the final result. 

Critical behavior corresponds to a nonanalytic dependence of the free energy in the 
coupling constants in V. For fixed couplings, this can be expressed as nonanalyticity in 
M//3 and P//3. Such behavior arises when the eigenvalue density acquires zeroes at the 
boundaries beyond those evident in (||). This saddle-point discussion is connected to the 
planar limit description (and eq. (||), in particular) by the relation | J\ — 4ABu(A)u(B) . 

We are now ready to calculate the free energy for the simplest non-trivial potential: 



Using eq.'s 



V(y) = —y + —y ■ 

and fli~0|) , the eigenvalue density is given by 



2ixM 



a 1 + ^(A + B] 



+ a 2 



(B-z)(z-A) 



(12) 



(13) 



We evaluate (|TT|), using ([IB]) and setting zq = (A + B)/2, to obtain 



P 2 



(B-A)< 
32 



ai 3 + ^ ai a 2 (A + B) + ^|(9/4 2 + 9B 2 + UAB) 



+ 8(3 
M 
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a x . A n . 2M P 
-(A + B) + — + - 



3b 2 + 3a 2 + 2ab - 4{a + b)vab 

b + Va 



P 



2M P 



log 



J, 



' B-A 
log — - — + 



The endpoints A and B are determined by 

a i 



2M P 
P 



log AB . 



a 2 



^(A + B) + ^(3A 2 ^ 

VabU + ^a + b] 



3B 2 + 2AB) 



(14) 

(15) 
(16) 



Given these expressions, it is not at all apparent which matrix potentials will yield nonan- 
alytic behavior in the free energy. Our approach is to determine the critical point from the 
planar string equations (|5]), and then examine the free energy (14) for nonanalytic behavior 
using the relations [pj] 



A=(V9-^fi) 2 B = (Ve- 

For the potential (0), the vanishing of the jacobian determinant becomes 



\J\ 



Setting g c 



a{ + 4aia 2 (# c + <f> c ) + Aa z 2 (6 z c + 6 c <p c + <p z c ) = 
6 C and C = y 2 with < y < 1, eq.'s (||[18D yield 

1 + y + y 2 1 2/ 3 (2 + 2/) 



a i 



2y + l 



(>2 



2y + l 



Qc 



+ 1; 



(17) 



;is) 



(19) 



One easily confirms that for these couplings the eigenvalue density p(z) has an extra zero 
at z = B = {l + yf. 

To study the singularity, one expands (|5|) around these critical values 



Aq 



2y + l 

2y 



(A0 - yA6) 



2y+l 



(A0 - yA6) 



2y 
l 



■(A9 2 + 2A9A<P) 



2y+l 



(A0 2 + 2A9A<f)) 



where A is used to denote Ax = x—x c . One observes that both Ag and Ag are proportional 
to A<j) = A0 — yA6, and that Aq = Aq + yAg has no linear variation. Further for generic 
rectangular matrices (i.e., < y < 1), it is convenient to define 



9 = 9- 



2y+l 



3^/(2/ + 1 

With these choices, the planar equations become 



9 = 6 + 



22/+ 1 
32/(2/+ 1)' 



Aq 



22/ + 1 
1 



A0 + 



22/ 



2(2/ 2 + 2/ + l) 
3t/(2y+l)(i/+l) 

y 2 + y + 1 



A9Acj) + 



(y-l)(y + 2) aX2 



9y 2 (y + l) 2 



-Acf 1 



3y{y + l)A9 



32/(2/ + 1] 



A0 2 



(20) 



The quadratic form of the singularity in (|20|) suggests Ag w — 5 2 t, where 5 = /3~ 2//5 
just as for quadratic singularities in double-scaling[l-3].Then eq. (|20|) yields Ag = —5s, 
AO = — J2 P 3 p fp an d Acf) = — J2 P 8 p h p , where f p and h p are solved for in terms of s and t. 
Upon inserting these scalings in eq. (|14]) , the planar free energy becomes 



122/(2/ + 1) 
5(22/ + l) 2 



(y 2 + y + l)(2y + l) 2 2 



36y 2 (2/+ l) 2 



s A + 



22/ + 1 
32/(2/ +1) 



-i 5/2 



(21) 
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This expression differs from the planar contribution to free energy presented for this model 
in ref. 0, only by rescalings of s and t, which were made there (and also in the next section) 
to simplify the results. Without such rescalings, both expressions would be identical. 
This demonstrates that the recursion relations successfully determine the correct critical 
behavior for the matrix model. Note that non-universal contributions to F were dropped 
from In fact, some of these terms are actually divergent, being proportional to /3 Q//5 
with a > 0, but they are analytic in both s and t. Similar divergent analytic terms 
appear in the scaling analysis of the saddle point approximation for the hermitian matrix 



model|]14[. A finite term analytic in s was also dropped. In general, one might expect 
additional finite terms which are nonanalytic in s0, but they do not arise for this critical 
point. 



3 Triple-scaling limit 



In this section, we reexamine the triple-scaling anaylsis of ref. j7|. We formulate the 
potential independent relations in terms of flow equations, and this allows us to determine 
the general string equation for this class of multicritical points. 

In these models since there are two parameters N and M which diverge separately as 
(3 — > oo, it is natural to begin with a general scaling ansatz for two linearly independent 
combinations of P and M 



AM + CP 

P 



AM + CP 

P 



p v ~H 



BM + DP 

P 



BM + DP 

P 



P 



0-1. 



In the scaling limit then, finite differences become derivatives: Dm = —P u Ad t — j3~ u Bd s , 
and dp = —f3~~ v Cd t — /3~ u Dd s . The planar limit analysis suggests v = 3v. As in the 
previous section, we set 9 C — 1 and <f) c = y 2 , and define 5 = (3~ v '. We then introduce 



e 



(p+i) 

M+n 



exp 



-8{An + Cl)d t - 5 3 (Bn + Dl)d s ] £ 5 q h q (t, s) 



q=2 



dP+l) 
} M+n 



y -exp 



-S(An + Cl)d t - 5 3 (Bn + Dl)d 8 \ ^ 6 q k q (t, s) (22) 

q=2 



First we determine the scaling limit of eq. (0), which are independent of the matrix 
potential. Inserting (|22|), leads to k 2 (t,s) = y(h 2 (t,s) + g 2 (s)) and a partial differential 
equation 



d,h 2 



C 2 y 



AT y 2 



C 2 y 



6 



y 



c 2 



2Ty 



: g 2 d t h 2 - -d s g 2 



(23) 



where T = AD — BC. It was noted in ref. 0] that neglecting the terms involving g 2 , eq. ( p3|) 
is the KdV equation with s and t playing the roles of the time and space coordinates, 
respectively. Presently, we also point out that eq. (|2~3D has the form of a flow equation. 
This connection is made clear by first introducing h(t,s) = h 2 (t,s) + g 2 (s)/2, and then 
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using (some of) the freedom to rescale s — > as, t —> f3t, h —> •yh, and g 2 — > V92- Thus 
eq. (p3|) can be expressed as 

d s h = d t (2g 2 (s)R 1 [h}+R 2 [h}) (24) 

where Ri[h] = —h/A and R 2 [h] = (3h 2 — d 2 h)/16 are the first and second Gel'fand-Dikii 
differential polynomials ||15| . 



Eq. ( p^) may be compared to the flow equations which are usually discussed in the 
context of matrix models and the KdV hierarchy|8|]. Given the string equation for the 
general model interpolating between multicritical points 

00 1 

t = J2( k +n)»kRk[h] , (25) 

fc=0 1 

the generalized KdV equations arise as flow equations for h({nk},t) 

dh 

d t R k+1 [h] . (26) 



If the first two couplings in the string equation (|25| ) were correlated as fix = s and = 
2g 2 (s), eq. ( p6|) would lead to a flow equation of the form given in eq. (|24"D. We reiterate 
that these results follow without restricting the matrix potential in any way. The only 
assumption is to fix the ratio of the scaling exponents, vjv = 3 (and q > 2). We have not 
found any other choices of this ratio which lead to interesting results. 

Next we would like to examine the string equations which arise from the potential 
dependent relations (||). For a quadratic potential (|TJ) with the critical values ( fL9l) and 
v = 1/5, one finds g 2 = s/2 and 

t = S ^R [h] + l 8 Ih[h] - I R 2 \h] (27) 

where Ro[h] = 1/2 and we have used the remaining freedom in rescaling variables to simplify 
these results. Comparing eq.'s ( p^) and (p5|) , one sees that /i and \i\ have precisely the 
dependence on s to be compatible with the flow equation fl24l) . 

Ref. also presents results for the critically tuned quartic potential which produces a 
third order multicritical point. Recasting those results in the form of eq. (^) by a rescaling, 
one finds that g 2 = (s/6) 1 / 2 and 

3 / \ 4 

t = 4 (ty R [h] + Is R x [h] - y f^J R 3 [h) (28) 

where R3IJ1} = — (10h 3 — lOhdfh — 5(d t h) 2 + dfh)/QA. We have also carried out the triple 
scaling analysis for the k = 4 multicritical point, which is produced by tuning a sextic 

5 This string equation {^) is the same as that presented in ref. but with slightly different rescalings, 
and a distinct combination of /12 and 172 for the scaling function. 
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potential. The results are similar to those above with g 2 oc s 1 ^ 3 and in eq. fl2"8p, s 3 ^ 2 — *> s 4//3 
and Rs[h] — > Ri[h], as well as changing numerical factors. For higher order critical points, 
we found it impractical to explicitly carry out the complete triple-scaling analysis, but the 
form of the general string equation is clear from the above examples. It takes the form 
of eq. (|25|) with the only nonvanishing coefficients being /i , Hi and /i^. The values of fi 
and Hi are fixed by the analysis of the potential independent equations (P^D- So it remains 
only to fix g 2 and The full triple-scaling analysis is unnecessary to determine these 
coefficients, but rather they can be extracted from the planar limit alone. At the (/c+l)'th 
multicritical point, we find 



92 



(k - l)\k\ 
2{2k- 1)! 



l/k 



and the string equation becomes 



^\Ro[h] + \s Rl [h] 



J/-1. 



2k 



Rk+i[h] 



(29) 



The string equation is that of the conventional (fc+l)'th order multicritical point but 
perturbed by the operators Ri and R . The new scaling parameter s governs the strength 
of the perturbations. 



4 Parametric solutions 

A seperate class of novel critical points is governed by coupled differential and finite dif- 
ference equations. These occur in the special limits where only a single large N parameter 
diverges. At present, we have no further remarks on the case of the vector models. We 
wish to point out the connection of the almost square matrix models, in which P — N — M 
remains a finite parameter, to the open-closed string equations proposed by ref. [§]. The 
simplest critical point was discussed in ref. J7| where the scaling function satisfies the 
Painleve II equation with a constant 

ld^-kr + skp=p + i 

(k[ and s will be defined below). Comparing to ref. suggests that P plays the role of 
an open string coupling constant. In these parametric models, the potential independent 
recursion relations give rise to a finite difference relation 

k[ P)2 - d s k[ P) = k[ P+1)2 + d s k{ P+1) ■ (30) 

The same relation was also developed for scaling functions satisfying a string equation 
in the mKdV hierarchy with coupling constants which differ by one|9|]. To confirm the 
connection of these parametric critical points to open-closed strings, we consider the next 
multicritical point. 
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This multicritical point can be produced by tuning a cubic potential, V(y) = + 
^fy 2 + ifl/ 3 - Square matrices correspond to g c = q c = 1, where the latter is a choice of 
normalization as is 9 C = 1. These multicritical points are related to approach of zeroes in 
the eigenvalue density to the origin. So first we choose <j> c = 1, which by eq. (|TT| ) yields 
A = and B = 4. Now eq.'s (||) and (H) fix ai = —2 + 2a 3 and = 1 — 4a3. Then, the 
eigenvalue density is given by 



Choosing 03 = 1/2 yields the form p(z) oc z 3 ^ 2 at the origin. 

For a scaling solution, we define 5 = and we use the following scaling ansatz: 
f = 1 - S^s and 

p=2 

Oih n -^ +N = -2e W (-n5d s )J25^(s) (31) 

P =i 

Derivatives do not replace finite differences in I because P does not scale. Given eq. ( |31~D 
without any further assumptions, the potential independent relations (0) yield eq. ( |3~0l) 
which relates the scaling function at different values of P. The identical equation is found 
in ref. [f| to relate solutions of the unitary matrix models with different numbers of flavours 
of 'quarks' U^]. The scaling of the potential dependent relations (|]) yields v — 1/5, 7 = 4, 



a 1 (P) 1 (P)2 

4 P) (*) = ^ (32) 



(P) 

and the scaling function k\ satisfies 

sk^ - M P)5 + \k[ P \d s kr? + \k^dlkP - \dtkT ] = P+\ (33) 

where we have shifted ki — ► aA^ and s — > s/ct with a = [8/3] 1 / 5 . This is the second 
equation in the mKdV hierarchy, where again P appears as the open string coupling in the 
constant term. 



Using eq. ([H]), one obtains the free energy 



OlF^ = -i(Af )2 + d s k[ P) ) (34) 



Here, the specific heat is proportional to the Miura map of the scaling function. Defining 
11 = k[ P ^ 2 + d s k[ P \ using eq. ( [331) may express the string equation in terms of the specific 
heat 

uR 2 [u] - ]-R[u\R"[u\ + ^(R'[u]) 2 = P 2 (35) 



9 



where R = R^iu] — s. This is the original form of the scaling equations proposed for 
open-closed strings with the open string coupling equal to P|| [L7|]. For P = 0, the 
perturbative solution of eq. ( |35D is simply R[u] = 0, which in the present case is the 
Painleve I equation. As well though, eq. provides a nonperturbative solution for two- 



dimensional gravity [17]]. The rectangular matrix models provide an explicit matrix model 
realization of these scaling equations extended to a nonvanishing open string coupling 
constant. 



5 Discussion 

Rectangular matrix models display a rich variety of multicritical behaviors. Underlying this 
diversity is the fact that these models have two independent large N parameters, which 
leads to multicritical behavior governed by two-dimensional singularities. 

Based on the analysis of the recursion relations alone, the connection of various critical 
points to singular behavior of the actual matrix model is not always clear. For the triple- 
scaling points, we have verified the validity of critical behavior by recovering the same 
nonanalytic planar contribution from a saddle point analysis. The multicritical points for 
the triple-scaling limit are governed by string equations which are simply expressed in terms 
of the Gel'fand-Dikii differential polynomials||15||. The usual fc'th order string equation in 



the KdV hierarchy is perturbed by the R± and Rq operators. The corresponding coupling 
/ii and /io are functions of the new scaling parameter, and a flow equation expressing this 
dependence naturally arises. The effect of the Rq can be absorbed in a renormalization of 
the cosmological constant to r = t + Xs^^" 1 ', where X is some numerical constant. For 
small s, as well as the usual genus expansion in r -( 2fc + 1 )A' ; there is also an expansion in 
s j T (k-i)/k a £ eacn g enus due to the perturbation of R\ in the string equation. In ref. 0, 
it was conjectured that the new dynamics uncovered by the triple scaling analysis should 
be related to a gas of punctures arising from the distinction between M and N loops in 
the surface interpretation In the present context then, it appears that tuning for triple 
scaling involves tuning these punctures to behave as a linear combination of the Rq and 
Ri operators. One might expect that a more subtle scaling would lead to more complex 
perturbations by linear combinations involving higher order Gel'fand-Dikii polynomials, 
but as yet we have been unsuccessful in producing such a tuning. 

The present reformulation of the triple scaling results provides some insight into the 
interpretation of the dual expansion for large \s\, which was noted in ref. J7|. The Ri 
perturbation expansion in \s\/r^ k ~ 1 ^ k is expected to have a finite radius of convergence, 
and beyond that point the solution should be expanded in terms of r/|s| fe ^ fc_1 ^.|18] In 



this domain, the model is expected to be in the neighbourhood of the "topological" model 
governed by Ri[h}. (This crossover behavior can be explicitly seen for the k = 2 and 3 
critical points, at least in the planar limit where the string equations are easily solved.) 
For this interpretation to be applicable, one should consider the real root of the planar 
equation which vanishes as \s\ — > oo. This leads to the expected expansion, h « -§7 + .. •• 
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Note that in this phase, the full expansion is not entirely a perturbation expansion in 
T l^ s \k/(k-i) which j s planar contribution of Rk[h], but rather the derivative terms in 
Rk[h] produce higher genus contributions as well. This is not the complete story though. 
For large \s\, the planar equation may have another real root which yields an expansion 
h oc s 1 /^ -1 ) + . . . (here if (k-1) is even, one requires s > 0, and then there is a second root 
at -s 1 /^" 1 )). These were the dual expansions presented in ref. J7|. In these cases, both Ri 
and Rk are equally important in determining the behavior of the perturbative expansion. 
The interpretation of these extra strongly coupled phases remains unclear. Similar strongly 
coupled domains may also separate the R\ and Rk phases, if there is a gap between the 
regions of convergence of the expansions in \s\/r^ k ~ l ^ k and r/\s\ k ^ k ^ 1 \ 

Finally, we have examined the parametric solutions for almost-square matrices. By 
examining the k = 2 critical point, we confirmed that these solutions provide an explicit 
matrix model realization of the scaling equations for open-closed strings 0, with P = N—M 
playing the role of the open string coupling constant. This is a natural position for P to 
appear in, since in the limit of almost-square matrices, the surface interpretation associates 
P with introducing punctures or boundaries^]. The scaling equations, which naturally 
arise with our scaling ansatz fl3~T|), are those of the mKdV hierarchy but with a constant 
P + |. The same equations naturally arise in the double-scaling limit of unitary matrices 
coupled to C = P + \ flavors of quarks |16|. We emphasize though that the physics of 



these models is not identical. In the present case, the free energy is simply related to the 
Miura map of the scaling function ([34]). For the unitary matrix models, one has simply 
Q 2 s P^ p ) = _I&( S ) 2 ; where k(s) is the corresponding scaling function[16|. Thus even though 



one has a simple map between solutions in the KdV and mKdV hierarchies^], the physics 
of unitary and rectangular matrix models remains distinct. 
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